Topology and energy transport in networks of interacting photosynthetic complexes 
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We address the role of topology in the energy transport process that occurs in networks of photo- 
synthetic complexes. We take inspiration from light harvesting networks present in purple bacteria 
and simulate an incoherent dissipative energy transport process on more general and abstract net- 
works, considering both regular structures (Cayley trees and hyperbranched fractals) and randomly- 
generated ones. We focus on the the two primary light harvesting complexes of purple bacteria, i.e., 
the LHl and LH2, and we use network-theoretical centrality measures in order to select different 
LHl arrangements. We show that different choices cause significant differences in the transport 
efflciencies, and that for regular networks centrality measures allow to identify arrangements that 
ensure transport efficiencies which are better than those obtained with a random disposition of the 
complexes. The optimal arrangements strongly depend on the dissipative nature of the dynamics 
and on the topological properties of the networks considered, and depending on the latter they are 
achieved by using global vs. local centrality measures. For randomly-generated networks a random 
arrangement of the complexes already provides efHcient transport, and this suggests the process is 
strong with respect to limited amount of control in the structure design and to the disorder inherent 
in the construction of randomly-assembled structures. Finally, we compare the networks considered 
with the real biological networks and find that the latter have in general better performances, due 
to their higher connectivity, but the former with optimal arrangements can mimic the real networks' 
behaviour for a specific range of transport parameters. These results show that the use of network- 
theoretical concepts can be crucial for the characterization and design of efficient artificial energy 
transport networks. 

PACS numbers: 87.15.hj, 82.39. Rt, 89.75.Hc, 71.35.-y, 05.40.Fb, 82.37.Vb, 87.15.A- 



I. INTRODUCTION 

Although research on energy transfer in photosyn- 
thesis has a very long history the refinement of 
experimental techniques in the last decade has dra- 
matically enlarged the range of possible observations 
at the molecular level 0, thus boosting new interest 
in the field. Recent progress in this context includes 
new experimental results aimed at characterizing on 
one hand the structure of light harvesting systems in 
other biological organisms |3| and on the other hand 
the presence and relevance of quantum effects in the 
energy transport process [ll|, [13 • These results are 
of fundamental interest since they allow to shed light 
onto the mechanisms at the basis of energy transfer 
processes, that owing to the long course of natural 
selection are likely to be optimally efficient. Therefore, 
a more complete understanding of their features has 
also a potential technological impact in providing useful 
benchmarks on how to engineer artificial light-harvesting 
systems. 

Purple bacteria are among the most important or- 
ganisms whose photosynthetic apparatus is currently 
studied. The basic actors in the energy transfer process 
within these membranes are two kinds of molecular 
photosynthetic complexes^ called the LHl and the LH2. 
The latter play the role of antennas, capturing the 
incoming photons and tunneling the resulting excitons 
to the former, that also act as antennas and furthermore 
contain the reaction centers (RCs) where charge separa- 



tion is eventually induced [g. 

The global dynamics in these models is highly dependent 
on the membrane architecture. In particular, recent 
studies have revealed the peculiar capability of these 
bacteria to adapt the structure of their photosynthetic 
membranes to the illumination conditions characterizing 
their growth process Q. While theoretical models of 
the transport process taking place in purple bacteria 
membranes have already been proposed [8l-[l(| and they 
are based on classical, Markovian master equations, so 
far only some basic model architectures were considered 
and a comprehensive study of the relevance of the 
membrane topology for efficient energy transfer has 
yet to be developed. This may be very important 
for applications: Current techniques like nanoimprint 
lithography already allow to form regular patterns of 
optically active LH2 complexes and future develop- 
ments in this direction might allow to arrange also other 
complexes, according to more complex designs. 
In this context, a basic and general question is the 
following: how is the transport process affected by the 
topological properties of the underlying structures? 
Are there design principles allowing to select the most 
efficient topologies? In the present work we propose 
a framework in order to address this questions. We 
take inspiration from the real biological membranes and 
simulate the energy transport process 0, on more 
general networks in order to relate the dynamics to 
different topological properties. 

Indeed, light-harvesting membranes can be regarded 
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as networks [2T| where the different photosynthetic 
complexes, such as LH1/LH2, play the role of nodes, 
and links represent the physical interactions (exciton 
hopping) between neighbouring complexes. One can 
thus vary the topological properties of the system by 
choosing different network architectures and by choosing 
different arrangements of the complexes on them. As for 
the architectures, we will focus on well-known regular 
structures like Cayley Trees and Regular Hyperbranched 
Fractals [l^ [T7l - [20| as well as randomly-generated ones 
(representing, respectively, presence or lack of control 
over the networks assemblage). Once a specific archi- 
tecture is selected, one has to address the problem of 
selecting specific arrangements of LH1/LH2 complexes 
on it. In this work we propose to use some of the 
most common network theoretical tools used in the 
description of dynamical processes on complex networks: 
centrality measures. The latter allow to rank the nodes 
of a network on the basis of their centrality with respect 
to it and we will therefore use them as guiding principles 
for the choice of different arrangements. 

The results of our analysis show that the topology has 
a remarkable influence on the dynamics. In particular, 
for the selected structures, there is a high sensitivity of 
the efficiency of the energy transport with respect to the 
choice of the arrangement of photosynthetic complexes. 
The use of centrality measures, in particular for regu- 
lar structures, allows to highlight specific arrangements 
that ensure a better efficiency with respect to a random 
arrangement; we benchmark these results by comparing 
them with those obtainable by using an optimization al- 
gorithm and with those pertaining to the real biological 
networks. As for random structures, the efficiency corre- 
sponding to a random arrangement is comparable with 
that achievable with the use of centrality measures; this 
highlights the robustness of the dynamical process with 
respect to the lack of control over LII1/LII2 arrangement 
in randomly-assembled structures. 

This work is structured as follows: In Sec. |TT] we 
introduce the basic notions on purple bacteria's light- 
harvesting membranes and the energy transport model 
used in the paper. In Sec. IIIII we introduce the types of 
networks and the main network-theoretical tools used in 
this work. In Sec. IIVI we discuss our main results on the 
relation between network topology and efficiency. Sec. |V] 
closes the paper with some final remarks. 



II. PHOTOSYNTHETIC ENERGY TRANSFER 
IN PURPLE BACTERIA 

A. Structure of the membranes 

The starting point of our analysis is the energy trans- 
port process that occurs in biological organisms such as 
purple bacteria. By virtue of techniques like X-ray crys- 
tallography and atomic force microscopy, the structure 




FIG. 1: Schematic depiction of basic photosynthetic 
complexes of purple bacteria LII2 and LHl; they both 

play the role of antennas, capturing the incoming 
photons. The LHls also contain the reaction centers 
(RC) where charge separation eventually occurs. 



of photosynthetic membranes of purple bacteria has been 
described with high precision Q ; moreover the main dy- 
namical timescales governing the exciton dynamics have 
been measured or calculated 0, [H. The energy trans- 
port process is based on two different kinds of pigment- 
protein photosynthetic complexes called LHl and LH2 
which grow in the membranes. The LH2 play the role of 
antennas, absorbing incoming photons and transferring 
the resulting electronic exciton to the LHl. The LHl 
complexes are also effective in absorbing light, at a dif- 
ferent wavelength, and they contain the reaction centers 
(RCs) where the exciton is absorbed by a special pair of 
clorophylls, inducing ionization of the pair - charge sepa- 
ration - and triggeringa chemical reaction (the reduction 
of quinone to quinol [8[). 

LH2 complexes are always more abundant than LHl 
ones. However, the stoichiometry s = N2/N1, i.e., the 
ratio between the number of LH2 {N2) and the number 
of LHl (A^i), varies depending on environmental con- 
ditions. The membranes have a remarkable ability to 
adapt to the intensity of the illumination during their 
growth: Bacteria grown under low light intensity (LLI) 
conditions have a stoichiometry s^li « 7 — 9, while those 
grown under high light intensity (HLI) conditions have a 
stoichiometry Shli ~ 3.5 — 5. As for photon capture the 
membranes' total cross section is cr = criNi-\-a2N2'- since 
0-2 = 116A^ > CTi = 67.29A^, the higher the stoichiomet- 
ric ratio, the higher the total cross section. The rate 
of photon capture R is proportional to cr: R = la/hv, 
where / is the incoming light intensity at the relevant 
wavelength. As a consequence, under a given illumi- 
nation instcnsity /, LLI membranes are able to capture 
more photons than HLI ones, which enables them to sur- 
vive in scarcer illumination. 
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B. The model of energy transfer 

Excitations arc initially created on one complex via 
photon absorption and can hop to neighboring complexes 
through a Forstcr resonance mechanism Q , until they ei- 
ther dissipate (with typical dissipation time tdiss) or lead 
to charge separation in a RC (with typical time tcs)- A 
schematic depiction of the process and the basic com- 
ponents is given in Fig. [TJ A crucial feature of the en- 
ergy transport process is the fact that when an exciton 
leads to charge separation within a RC, the latter fea- 
tures a "busy" time interval, during which quinol is pro- 
duced, removed and then a new quinone becomes avail- 
able. Therefore the RC is closed for a time interval tbiock , 
called recycling time, during which it cannot exploit in- 
coming excitons. Furthermore, the RC in the closed state 
becomes a quencher for excitons, i.e., increases its dissi- 
pation rate (t^J^ ^ tdiss)- )■ The timescales for energy 
capture, exciton lifetime and RC reopening are roughly 

Ims respectively. The recycling time tbiock can differ 
significantly depending on the specific kind of bacterion 
analyzed (1 — 30ms). As a consequence of RC closure, 
at any time only some of the RCs are available for pro- 
ducing charge separations. The number of available RCs 
decreases with increasing tuock - i-C, as the RCs remain 
closed for a longer time - and the overall transport pro- 
cess passes from an active regime in which all RCs are 
open to a saturated regime in which all RCs are closed. 
In Table HI and HIl we report the values of the relevant dy- 
namicalparameters of the LLI- and HLI- adapted mem- 
branes 

The problem of modeling energy transfer processes 
in biological and artificial systems has been thoroughly 
studied in the past decades, and models of the exciton 
transport dynamics in light-harvesting membranes have 
been developed in order to take into account both inco- 
herent and coherent phenomena [l^l- Here we focus on 
a simple classical, Markovian model developed in 
to model energy transfer in purple bacteria membranes. 
Since exciton transfer between complexes in such mem- 
branse arises through the Coulomb interaction on the 
ps time-scale, while vibrational dephasing destroys co- 
herences within a few hundred /s, coherent effects in 
the energy transfer between complexes are expected to 
be very weak and are therefore neglected in these mod- 
els [HI 

Since tdiss ^ R^^, simultaneous occurrence of two exci- 
tons in the membrane is very unlikely, therefore the sys- 
tem can be modeled as if a single exciton were present at 
each time t. The dynamics is described by a Markovian 
Master equation (see appendix B for details) and thus 
can be numerically simulated by standard random walk 
methods [§, [l^l ■ In the following, we will always base our 
analyses on random walk simulations. 
Excitons are created at random times {ti} determined in 
advance by using a Poissonian distribution. At each time 
ti a single exciton is randomly created in a LHl (LH2) 



TABLE I: Relevant scale parameters for HLI and LLI 
networks: stoichiometric ratio, absorption cross-section, 
and dynamic timescales 
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TABLE II: Typical transfer times between different 
photosynthetic complexes 
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exciton then follows a random walk, and the probability 
of jumping from a site j in a given time step St is given 
by Pj„™p(i5t) = KjSt with 
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exciton transfer rates Wij are zero if the nodes are not 
linked, else they are taken as the inverse of the typical 
exciton transfer times tx^Y- The values of all timescales 
are the same as in Tables U and [Hi The excitons remain 
on the same site or jump to a neighboring one according 
to different probabilities, until they dissipate or lead to 
charge separation in a RC. When charge separation oc- 
curs in a RC, the latter remains closed for time tbiock in 
the simulation and its dissipation rate is increased. Our 
simulations cover a long time of ^ Is, so that the histo- 
ries of Nabs = 10^ excitons are simulated [iGj . 
The optimality and robustness of energy transfer are 
characterized by performance measures: The most rel- 
evant is the efficiency of the transport process 77, i.e, the 
probability that an exciton lead to a charge separation 
- while 1 — 77 is the probability that it dissipate. In the 
random walk, we evaluate the ratio r] = Ncs/Nabs, where 
Ncs is the number of photons that are used for charge 
separation in RCs and Nabs the total number of photons 
absorbed by the network. 



III. BIOLOGICALLY-INSPIRED NETWORKS 
AND NETWORK-THEORETICAL TOOLS 

In order to study the role of topology for the biologi- 
cally inspired transport process described in the previous 
section, we will consider different kinds of regular and 
randomly generated networks. The structures that we 
are going to use are based on a pair {F,^} composed 
by a graph F = {E,V), with set of vertices V and 
edges E, together with an arrangement A = {Vi,V2) 
i.e.; a subdivision of the vertices in two disjoint subsets 
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Vi,V2 C V that identify the nodes that are occupied by 
LHl and LH2 respectively. The complete structure is 
obtained by adding to F one node for each element of Vi , 
representing a RC directly linked only to the respective 
LHl complex. We have l^il = iVi and |]^| = A^2, where 
the number of the different complexes is determined 
by the stoichiometry s. We will investigate two differ- 
ent settings, representing specific high light intensity 
(HLI) and low light intensity (LLI) conditions. The 
stoichiometric ratios are those found in real networks: 
N2/N1 = 4.64 for for HLI, N2/N1 = 7.04 for LLL 

In general, the transport process properties will be 
determined both by the topological features of the 
graphs F and by choice of the arrangement A, i.e, 
how the LH1/LH2 complexes arc placed on F. The 
different arrangements A will be chosen following criteria 
based on tools developed in complex networks theory 
(centrality measures). The selected structures will be 
compared with the real networks describing purple 
bacteria membranes, from the point of view of energy 
transfer properties. 



A. Networks 

The regular networks we will focus on are well-known 
regular structures which have been considered [T7l - [l9| 
in models of energy transport: i) C'ayley trees (CT) ii) 
regular hyperbranched fractals (RHP). 
The CT are constructed as follows. The n-th generation 
d-Cayley tree is a tree of n levels in which all vertices 
on the interior have degree d, while vertices on the out- 
ermost layer have degree 1. In our simulations we use 
d = 4, n = 4 (see fig. ^ and a total number of nodes is 
161. 

The RHF instead are constructed as follows. A first 
generation (RHF) (as shown in fig. [2]) of functionality 
/ is a star graph consisting of a central vertex connected 
through / edges to / surface vertices. To construct a sec- 
ond generation RHF, / copies of the first generation RHF 
arc connected to the core first generation RHF through a 
single leaf-leaf edge. This procedure is repeated n times 
for an n-th generation RHF. In our simulations we use 
/ = 5, n = 3 and the total number of nodes is 216. Both 
CT and RHF depicted in fig. H 

In addition to regular structures, we also investigate 
randomly-generated ones. This allows to compare the 
effiency of transport achieved in two very different cases: 
the case where one has global control over the network 
topology F and the case where one has no control on 
it - this might happen, for instance, if the network is 
randomly assembled. A natural way to generate pla- 
nar structures that are possible realizations of light- 
harvesting networks is to consider for example Randomly- 
decimated hexagonal Networks (RN). These are obtained 
according to the following prescription: We consider a 
hexagonal lattice of Ntot = 256 sites and remove m 




FIG. 2: Regular graphs, generated according to the 
procedure in subsection IIIIAI (Top) Cayley tree with 
d ^ 4, n 4 (Bottom) RHF with / 5, n = 3. In both 
graphs, nodes of different colors belong to different 
generations: (yellow), 1 (red), 2 (green), 3 (blue), 4 
(black). 




FIG. 3: Random network generated according to the 
procedure in subsection IIIIAI 



randomly-chosen links such that the connectedness of the 
network is preserved. In our simulations m = 193, 321 so 
that the average degree becomes ^ 4 and ~ 3 respec- 
tively. For both values of m we generate an ensemble of 
M = 100 RNs (a particular element of the ensemble is 
depicted in fig. |3]). 

Finally, we consider the real, biological networks stud- 
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FIG. 4: Simplified depiction of a typical HLI 
photosynthetic membrane of purple bacteria 



ied by Fassioli et al. [0| , which are based on studies of the 
Rsp. photometricum membranes 01 . The global mem- 
brane structure is known in great detail, since atomic 
force microscopy has provided high-resolution images of 
the membranes of the purple bacterium Rsp. photo- 
metricum 0- The membranes have a planar structure 
where complexes are densely packed, and the cores are 
surrounded by LH2 complexes (5-7 on average) 0- A 
simplified depiction of a typical HLI membrane in terms 
of the constituent complexes is given in fig. UJ in the 
network model, links are present between neighbouring 
complexes. Notice that we consider a single realization of 
real network for HLI and LLI i.e., in both cases F and A 
as fixed and correspond to the networks in Q . Further- 
more, LLI and HLI membranes, considered as networks, 
have a different topology: since LHl have on average 
more links to neighbouring complexes than LH2, HLI 
structures have a higher average connectivity. The total 
number of nodes are 168 (HLI) and 193 (LLI). 

The parameters defining all the above networks have 
been chosen in order to have structures with approxi- 
mately the same size i.e., the same number of nodes \V\. 



B. Choosing different arrangements: Centrality 
measures 

Once a specific network F is chosen, another step is 
needed in order to completely specify the topology of 
the transport network. Indeed, different arrangements 
A of the LH1/LH2 complexes are possible and one has 
therefore to identify guiding principles allowing for the 
selection of different incquivalent arrangements. In this 
sense, a possibility is given by the study of the relevance 
of the nodes within a network. Great efforts in the de- 
velopment of the theory of complex networks have been 
devoted to this task and in this context several central- 
ity measures have been proposed (22j . The usefulness of 
a given measure depends strongly on the specific prob- 
lem at hand, since the latter determines which nodes are 



more relevant from case to case. 

Two common families of centrality measures are the set 
of degree-like and the set of betweenness-like central- 
ity measures. Degree-like measures are essentially based 
on counting the number of paths emanating from the 
node. Different measures are obtained considering differ- 
ent path lengths and different kinds of paths (edge- of 
node- disjoint, geodesic, etc.). The simplest of such mea- 
sures is the degree(DF,G) [2j]: The degree fc of a node 
is simply the number of nodes to which it is connected 
(i.e., the number of paths of length 1 emanating from the 
node). Betweenness-like measures instead are based on 
counting the number of paths that pass through a given 
node k. Again, different measures are obtained consid- 
ering different path lengths and different kinds of paths. 
The most common of such measures is the shortest path 
betweenness centrality (SPB) [g^j; here the basic objects 
are the shortest paths connecting two nodes i,j (a path is 
given by a sequence of edges connecting the two nodes) 
; in general there are many different paths of minimal 
length (sequences with minimal number of edges) be- 
tween two nodes. The SPB of a node k is defined as 
Bk ~ where gij is the total number of shortest 

paths from node i to node j, while gijk is the subset of 
such paths passing through node fc. While the degree is 
a local measure of node centrality, depending only on the 
structure of the network in the immediate neighbourhood 
of a given node, the shortest-path betweenness is a global 
measure of node centrality, since it takes into account the 
structure of the whole network. 

By means of the above measures one can identify dif- 
ferent incquivalent arrangements A in the following way. 
For each F we evaluate the DEG, SPB of all nodes (a full 
characterization of F^j^, Tjjhf F/^jv in terms of these 
measures is presented in the Appendix A). An arrange- 
ment A is identified by choosing the subset Vi as the set 
of A^i nodes satisfying one of the following criteria (s^ 

• nodes with maximal SPB, or DEG 

• nodes with minimal SPB, or DEG 

• nodes chosen at random 

While other choices are possible, these are the basic 
and natural ones that allow to test how the transport 
efficiency can be affected by different arrangements of 
the RCs on a given network. Like in the case of ran- 
dom topology, the random arrangement represents the 
case when no control over the LH1/LH2 arrangement is 
possible. 



IV. ENERGY TRANSPORT IN BIOLOGICALLY 
INSPIRED NETWORKS 

In the following we give a detailed description of the 
dynamics and then we discuss the relation between the 
dynamics and the topological features of the selected 
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structures. In order to describe the dynamics we intro- 
duce the foUowing functions of the recycHng time t^iock '■ 

• the efficiency rj = Ncs/Nabs] 

• the average fraction of closed RCs {nf^fj) = 

• the average lifetime (t)^^ of successful excitons, i.e., 
those that reach a RC and induce a charge separa- 
tion; 

• the average maximal distance {d)^^ of successful ex- 
citons, i.e., the distance, in steps, between the ini- 
tial site and the most distant site reached during 
the random walk); 

• the average exploration parameter (x)^^ of success- 
ful excitons, i.e., the fraction of network sites vis- 
ited. 

In particular (d)^^ and (x)^^ characterize the diffusion of 
excitons over the network. 

For each choice of T and A, and for any fixed value of 
tbiock we realize the random walk dynamics corresponding 
to Nabs = 10^ excitons. The above averages are taken 
over the Nabs histories (ry, (jt-ac)) or over the histories of 
the Ncs successful excitons ((t)^^ , (d)^^ , (x)^^). In case of 
random networks or arrangements the functionals are the 
results of a further average over the elements of a given 
ensemble of networks and/or arrangements. Finally we 
let tbiock varying over a large range of values (10~^ms — 
lO^ms). 

A. Cayley trees 

We now show how the introduced functionals allow 
to characterize the dynamics and its relation with the 
topology of LHl/RC arrangement. We shall first focus on 
CTs, and then analyse results for the remaining network 
topologies (RHF, RN). 

Different transport regimes. As a first general com- 
ment, we note that the dynamics is characterized by three 
main time scales, given by the absorption rate R~^, the 
typical dissipation time t^iss and the recycling time tbiock ■ 
In particular there are two distinct regimes: The active 
regime, characterized by tbiock < < tdiss, i-e., by 

a recycling time smaller than the absorption rate; and 
the saturated regime, characterized by R^^ < tbiock "C 
tdiss, i-e., by a recycling time greater than the absorption 
rate. The transition between the two regimes occurs at 

R^ ~ tbiock- 

In fig. [5] (upper panel) we show rj for CTs for all different 
LHl arrangements in LLI conditions. The behaviour for 
HLI networks (not shown) is analogous; since Rjj\j ~ 
10~^R2l^j the transition between the two regimes occurs 
in HLI networks earlier than in LLI ones. 

Each absorbed photon moves in a network which is 
characterized not only by the stoichiometric ratio, but 
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FIG. 5: (Top) Efficiency 77 of exciton transport for CT 
as a function of tbiock in LLI conditions, under different 
LHl arrangements. One can identify the active (shaded 

green) and saturated (shaded red) regimes. The 
transition occurs at tbiock ~ ~ 0.1 ms. (Bottom) ?/ 
and RC occupation n'j^^ for CT as a function of tbiock in 
LLI conditions, rj sinks and n'j^^j grows when tbiock is 
increased. For any fixed tbiock, efficient arrangements 
correspond to higher values of both rj and n'j^Q . 

also by the number of closed RCs. In fig. [5] (right panel) 
both 77 and n'^Q are plotted together for the different 
LHl /RC arrangements. On one side, for any given ar- 
rangement, grows with increasing tuock, while ry 
sinks. On the other side, if we compare different arrange- 
ments at a fixed value of tbiock, we notice that more ef- 
ficient arrangements have a larger fraction of closed RC: 
excitons are more likely to be absorbed, so that on aver- 
age more RCs are closed and consequently the ordering of 
the efficiency curves is the same as the RC-closure curves. 

As for the other functionals, (r)^^ , (d)^.^ , (x)^^ we note 
that they all display the same behaviour with tbiock- they 
remain essentially constant in the active regime and they 
grow throughout the transition region. Indeed when 
many RCs are closed, excitons are expected to travel fur- 
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ther away from their initial site and explore a larger frac- 
tion of nodes in order to find open RCs. A maximum is 
reached when nj^^ ^ 90%. Thereafter, unless one of the 
very few open RCs is found in the vicinity of the initial 
exciton site, the excitons are dissipated before reaching 
any RC and thus all quantities undergo a slight decrease 
after this threshold. 

The concentration of RCs is also at the basis of the differ- 
ence between LLI and HLI networks. Indeed, the values 
of (t)^^ , and (x)^^ are higher in the former case; 
due to lower stoichiometric ratio, the excitons have to 
explore a bigger fraction of the network; for instance, in 
the active photosynthesis regime one has that for LLI 
(r)^^ e [90ps,180ps] while for {tcs} £ [80ps,120ps] for 
HLI networks. 

The plots of (x)^^ in Fig. [H] (panels c and /) show that 
the absorbed excitons do not explore the whole network, 
but only a relatively small portion of it: only 5 — 10% of 
all sites are visited by successful excitons. Furthermore, 
the ordering of the (d)^^ and (tcs) curves is reversed with 
respect to the rj ones. We thus reach a first important 
conclusion: due to dissipation effects, excitons cannot ex- 
plore the whole network, and in efficient configurations 
excitons find open RCs in the vicinity of the initial site 
and reach it in the shortest time. Optimal arrangement 
criterion and topological properties of the network. In the 
case of the Cayley tree the LHl arrangement criterion 
giving the highest efficiency is the maximal SPB, and 
the corresponding subset Vi is composed by the nodes 
clustered around the root of the network. This fact can 
be explained by analyzing the topological properties of 
the network Tct and the dynamical behaviour. Indeed, 
in the active photosynthesis region the distance (d)^^ « 4 
is comparable with a key topological property of the net- 
work, namely, the average distance between two nodes of 
the network which is of 6.4 steps. This means that, al- 
though the exploration parameter is relatively small, the 
exciton is able to move across the network e.g., from the 
periphery to the center of the tree and vice- versa. There- 
fore, it can take advantage of the global character of the 
centrality of the nodes occupied by the LHl/RCs. This is 
confirmed by the positive value of the linear correlation 
coefficient between the SPB of each node and the aver- 
age time spent on it, which we find to be approximately 
0.7 in all regimes: excitons are likely to reach the central 
region and to spend more time on the Vi nodes. 
A general important feature highlighted by our analysis 
is that both for LLI and HLI structures the maximal SPB 
criterion allows for a significant improvement with respect 
to the random arrangement. This indicates that the cho- 
sen network theoretical tools represent a useful guiding 
principle for the identification of the most efficient con- 
figurations. 




(a) lifetime (t)^ 
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FIG. 6: Lifetimes (t)^^, average maximal distances 
(d)^^, exploration parameter (x)^^ of successful excitons 
as a function of the recycling time tuock for Cayley trees 
in HLI (left) and LLI (right) conditions, under different 
LHl arrangements: (black) random (red) max DEC 
(green) max SPB (blue) min DEC (yellow) min SPB. 
All quantities grow in the transition region, due to the 
necessity of longer paths in order to find open RCs, and 
have larger values for LLI networks that have less RCs. 



B. Regular Hyperbranched Fractals. 

In the case of RHF networks, the global behavior of 
all functionals with tuock is analogous to the one de- 
scribed for the CTs. However, some fundamental dif- 
ferences emerge. On one hand, the range of the various 
functionals is much higher for the RHF and this is a clear 
indication that these networks are very sensitive to the 
choice of the LHl arrangement. For example, in the ac- 
tive photosynthesis regime (t)^^ G [lOOps, 375ps] for LLI 
while for HLI (r)^^ e [70ps, 240ps]. 

On the other hand, the criterion that allows for a maxi- 
mal efficiency is in this case the maximal DEC, i.e., a local 
measure of centrality of the nodes, wheareas the arrange- 
ment determined by the maximal SPB criterion (Fig. [71 
left panel) has a considerably low efficiency (60%). This 
result can again be explained in terms of the topolog- 
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FIG. 7: (Top) efficiency curves for RHF at LLI under 
different LHl arrangements. Tlie arrangement based on 
global SPB is particularly inefficient. (Bottom) 

efficiency curves for RHF at HLI. Arrangements are 
based on max SPB{d), d = 2 — 6); efficiency grows as an 
increasingly local structure of the network is considered. 



ical features of the Trhf- Indeed, the values of (d)^^ 
are in general small if compared to the average distance 
between the nodes of the Thhf network which is 16.68 
steps (for example {d)^^ e [4, 6] for LLI case in the active 
photosynthesis region). Therefore, the dynamics can not 
exploit the global centrality of the nodes with max SPB 
i.e, those clustered around the root of the network; This 
fact is again confirmed by the linear correlation coeffi- 
cient between the betweenness of each node and the av- 
erage time spent on it which is now negative and equal 
to —0.5 in the case of max SPB arrangement. 
The previous result suggests a general consideration: for 
networks with high values of the average distance be- 
tween nodes and in presence of a dissipativc dynamics 
that does not allow for high values of (d)^^, optimal LHl 
arrangements can be found by considering local centrality 
measures, able to characterize the centrality of the vari- 
ous nodes with respect to the local surrounding network. 
In order to further test this fact we introduce a class of 
centrality measures which generalize the concept of local 
betweenness centrality already introduced in the context 



of search problems on complex networks [27|. For each 
node A: of a given network F, one can fix a distance d and 
consider the subnetwork composed by those nodes whose 
distance from k is < d. One then evaluates the short- 
est path betweenness of the given node in the selected 
subnetwork. We shall denote such "d-local" between- 
ness centrality measure as SPB{d) [1^. The latter coin- 
cides with the SPB for sufficiently high values of d). For 
each value of d the arrangement we consider corresponds 
to placing the LHl complexes on the nodes with higher 
values of SPB{d). Let us consider RHFs and compare 
the results obtained with the LHl arrangement based on 
global SPB and that based on SPB{d) for various val- 
ues of d. In Fig. [7K right panel) we clearly see, that while 
the global SPB induces a very inefficient arrangement, 
the local SPB(d) measures progressively induce more ef- 
ficient arrangements as d decreases and a increasingly 
local structure of the network is considered [l^. Even 
though this precise hierarchy is specific to the RHF, it is 
also a good qualitative example of how local, rather than 
global centrality measures may be more appropriate for 
dissipativc dynamics that take place on networks with 
large average distance between sites. 

We conclude our discussion by observing that also in 
the case of RHF structures, the use of criteria based on 
centrality measures allow to identify optimal arrange- 
ments whose efficiencies are significantly higher than 
those obtained by randomly placing the various com- 
plexes. 



C. Random networks 

We now pass to examine the random networks gener- 
ated according to the procedure described in Sec. IIIII 
Notice that in the following discussion and figures, we 
always cosider averaged quantities over an ensemble of 
M ~ 100 randomly generated networks. In Fig. [8] we 
show the efficiency curves for the LLI networks with 
TO — 193,321 links eliminated (average DEG= 3,4 re- 
spectively). 

The efficiency ?/, as well as the other functionals, glob- 
ally follow the same behavior of the previously exam- 
ined cases. The main difference between the two cases 
m = 193 and m = 321 is the global connectivity which is 
higher in the former case. This allows for higher values 
of {d)^g and {x)^^ and in general for slightly higher values 
of efficiency. 

In both cases, the criterion that gives the best efficiency 
is a local one, i.e., the maximal DEG. As for the rela- 
tion between topology and arrangement criteria, we first 
focus on the LLI case with to = 321, which are the net- 
works which show the maximal variation of the various 
functionals with respect to the change of the arrange- 
ment. Here the average maximal distance traveled in 
the saturation regime is (d)^^ € [5.8,6.5], which has to 
be compared with the average distance between pairs of 
sites which is « 11. This is consistent with the fact that 
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FIG. 8: efficiency for random networks at LLI under 
different LHl arrangements: (Left) m = 193 (Riglit) 
m = 321. A lower connectivity leads to a slightly lower 
efficiency. The best topological criterion (max DEG) 
has approximately the same efficiency of the random 
arrangement. 



the criterion giving the best efficiency is a local one, i.e., 
the maximal DEG. However, for the case m = 193, the 
values of {d)^^ are in general higher than the m = 321 
case while the average distance between pairs of sites is 
smaller (~ 9), which would suggest a global criterion like 
max SPB. Therefore we have an indication that in the 
case of randomly generated networks these two parame- 
ters are not sufficient to identify an optimal arrangement 
criterion. 

Moreover, at variance with the case of regular networks, 
for all values of tuock the efficiency corresponding to the 
best topological criterion (max DEG) is only slighty bet- 
ter than that obtained with a random arrangement of 
the LHls. This result is very relevant since it suggests 
that the dynamical process (originally modeled for pur- 
ple bacteria membranes) is a sense optimized to take 
into account the disorder inherent in the construction 
of randomly-assembled structures. In particular, this 
should hold for the real biological structures, but also 
for artificial structures where no control on the network 
structure F is possibile. Global optimization of the topol- 



ogy is indeed a rather difficult task (biologically and arti- 
ficially). Therefore, on one hand a dynamical process like 
energy transfer has to be strong with respect to limited 
amount of control in the structure (F, A) design, in order 
to be successful. On the other hand, if at least a local 
control in the topology of the network is possible (this 
should be in principle an easier task to accomplish), then 
the process could be devised in order to take advantage 
of the max DEG arrangement. 



D. Sensitivity of the efficiency to arrangement 
criteria in different networlc topologies 

As we have noticed, different LHl/RC arrangements 
can cause significant differences in efficiencies, lifetimes 
and all other functionals. In the following we shall focus 
on r], which is by far the most important one - being the 
main goal of biological and/or artificial optimization. In 
particular, the basic feature we want to study is the sen- 
sitivity of the efficiency of a given network with respect 
to the choice of the LHls/RCs arrangement. In order 
to quantify this sensitivity we introduce the functional 
Sr] = \r]M - Vml/VM where, for any fixed values of tuock, 
Vm iVm) is the maximal (minimal) efficiency correspond- 
ing to the optimal Am (worst. Am) arrangement. 
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FIG. 9: Sr] vs twocfe for the Cayley trees (left) and RHF 
(right). In LLI networks the LHl fraction is lower, so 
networks are more sensitive to LHl arrangement 

In fig. [S] the sensitivities of CTs and RHFs are shown 
as a function of tuock- A first important feature is that 
the LLI networks are in general more sensitive than the 
HLI ones. Thus, for higher values of the stoichiometric 
ratio, the definition of criteria for choosing RCs arrange- 
ments becomes crucial in order to have efficient energy 
transport. 

As for CTs the sensitivity has some initial value (5 — 9%) 
in the active regime, it grows to a maximum (w 12 — 14%) 
approximately when tuock = i ~ 2R~^ and then falls off 
to zero when tuock 3> R^^. It is therefore in the transi- 
tion regime that in general the energy transport is most 
affected by the choice of RCs arrangement. 

This feature is rooted in the ability of most efficient 
configurations to ensure high value of efficiency and at 
the same time a higher resilience to RCs closure. In 
order to clarify this point we focus on the absolute sen- 



10 



0.0-1 0.4 




1E-3 0.01 0.1 1 10 100 1000 10000 1E-4 1E-3 0.01 0.1 1 10 100 100010000 



FIG. 10: Relation between A77 and for Cayley 
trees at LLI: (left) drj/dn'j^Q (right) n'j^^ 



sitivity A77 ^ tjj^j — rj„i (here M = maximum SPB, m = 
minimum degree), which mirrors the behavior of ^77, and 
on its derivative with respect to the recycling time tuock ; 
in the following we define drj/dtuock = V- Iri the initial 
part of the transition region {tuock < t) we have that 
'Hm < "Hm < , i.e., the efficiency for the best arrange- 
ment Am decreases less than that of the least efficient 
one Am and therefore A?) > 0. The derivative of the ef- 
ficiency with respect to tuock can be decomposed in the 
product of two contributions f] = drj/dn'j^fj where 
is the number of RCs closed for a given arrangement. 
Since for all iwocfc one has > ^ > 0, the 

enhancement of the absolute sensitivity must be traced 
back to the variation of the efficiency with respect to 
the number of closed RCs; indeed for ttiock < t wc have 
drjm/dn'j^Q < dr]]\j /dn'j^Q < 0. This relation shows that, 
at least in the first part of the transition region (t corre- 
sponds to ~ 50%) the maximum SPB configuration 
is not only the most efficient but it is also more resilient 
to RCs closure. This feature no longer holds ior t > t 
and therefore when > 50% the sensitivity starts to 
decrease: > dijm/dn'^^ > drjM / dn'^^ and Afj < 0. 

As for the sensitivity of RHFs, we first notice that the 
initial sensitivities are very high: 21% in the HLI case and 
37% in the LLI case, while the enhancement of the sensi- 
tivity in the transition region is practically absent both in 
the HLI and LLI case. This is due to a lower resilience to 
RCs closure in the case of the best arrangement (maximal 
DEC) and also has an impact on the overall performance 
of RHFs. Indeed a direct comparison shows that the most 
efficient configuration for RHF is much less resilient than 
for CT, i.e., \drjM /dnj^(j\RHF < \drjM / dn^j^clcT and this 
causes the RHF to be the network with the worst perfor- 
mance in the transition region (see also fig. [T^ where ar- 
tificial networks' performances are compared with those 
of the real biological networks). 

We finally analyze the sensitivity of RNs. We have that 
Sr] for LLI varies between 5 — 10% (see Fig.[TT]), and its be- 
havior with tbiock is similar to the Cayley trees case. The 
parameters that influence the sensitivity are both the sto- 
ichiometry and the average connectivity determined by 
the number m of links eliminated from the original hexag- 




FIG. 11: sensitivity for random networks (Left) 
m ~ 193 (Right) m = 321. A lower connectivity 
enhances the sensitivity 



onal lattice. Indeed, a smaller connectivity (m ~ 321) 
makes the choice of the arrangement more crucial and 
therefore it enhances the sensitivity of the network, and 
this is true in particular in the transition region. 
To sum up, the high values of the sensitivity by all net- 
works analyzed indicate that the analysis of the topol- 
ogy of RCs arrangements becomes crucial in order to 
maximize the transport efficiency for such kind of net- 
works, in particular for LLI networks. The behavior of 
the sensitivity with respect to tuock is determined by the 
higher/lower resilience of efficient configurations to RC 
closure. 



E. Comparison with real networks and optimized 
arrangements 

In order to evaluate the optimality of the arrangements 
defined on the basis of centrality measures, we now com- 
pare the efficiency of all structures considered (CT, RHF 
and RNs) on one hand with those structures describing 
real, biological networks with a similar number of nodes 
and on the other hand with the results that one can ob- 
tain by searching for optimal arrangements by means of 
a numerical optimization method (annealing). 

To this aim in fig. ll2l and ll3l we plot Sx = {Vx—Vm) I Vm 
where x ={real, annealing} and tjm is the efficiency for 
the best arrangement of a given network. A first rele- 
vant message of these plots is that the real biological net- 
works are optimal both at LLI and HLI, in the sense that 
they have an equal or better performance than all other 
network topologies considered, and the main differences 
arise in the first part of the transition region. The opti- 
mality of real networks has to be related to their higher 
degree of connectivity with respect to the other networks 
analyzed. In particular they have a higher value of the 
average degree (< k >~ 5.75 for both HLI and LLI), 
and smaller values of the average distance between pairs 
of nodes (5.1 for HLI and 4.7 for LLI) and these features 
allow the excitons to explore greater parts of the network 
in search for an open RC; indeed, while for such networks 
< d >csG [5.5,7] is comparable or slightly higher than 
in the other cases , the exploration parameter reaches 
sensibly higher values: < Xcs >£ [0.16 — 0.29] for HLI 
and < Xcs >G [0.17 — 0.28] for LLI, where the minimum 
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corresponds to the active photosynthesis region and the 
maximum to the transition one. 

Despite their lower connectivity, however, there are pairs 
of artificial network topologies and RCs arrangements 
that allow for efficiencies that are comparable with those 
of the real networks at least in the active region; this 
in particular is true for Cayley trees for LLI conditions 
with LHl/RC disposed in nodes of maximal SPB and 
RHF for HLI conditions and LHl/RC disposed in nodes 
of maximal degree. As for the RNs, those that better 
approximate the efficiency of the real ones are those with 
a higher degree of connectivity (m = 193). 

Our analysis of artificial networks allows to identify the 
optimal pairs {F, A} that insure the realization of high ef- 
ficient energy transport in all regions: (Cayley trees, max 
SPB), (RHF, max DEC), (random networks, max DEC). 
In order to further test the optimality of such arrange- 
ments based on centrality measures, we compare their 
efficiencies with those obtained through an optimization 
method. For the latter, we have used an optimization al- 
gorithm based on simulated annealing [30| to search for 
the best A. The algorithm uses 1 — as a cost function, 
and evaluates it in the Master equation approach (see 
appendix B) in the limit of no RC closed {tuock ^ R~^)- 
Starting from a random arrangement, LHl and LH2 po- 
sitions are swapped till a minimum of the cost function is 
found. The optimal arrangement Aanneaiing is then used 
for all values of tuock- 

The comparison (sec fig. I13p shows that, in particular 
for regular topologies (CT for HLI and LLI, RHF for 
HLI) the best LHI/RC arrangements based on central- 
ity measures are optimal in that they allow to obtain, in 
all regions approximately the same efficiencies obtained 
by numerical optimization, while in the case of RNs the 
optimization procedure gives only slightly better results. 




FIG. 12: efficiency Srjreai vs. thiock for different 
networks, HLI (right) LLI (left). The real biological 
networks have a comparable (active region) or better 
(transition region) performance than all other network 
topologies considered. 



V. CONCLUSIONS 

In our work we have taken inspiration from biological 
light harvesting networks of purple bacteria and we have 
simulated an incoherent dissipative energy transport 




FIG. 13: efficiency 5r]anneaiing vs. tuock for different 
networks, HLI (right) LLI (left). For regular topologies 
the best LHI/RC arrangements based on centrality 
measures are optimal in that they allow to obtain 
approximately the same efficiencies obtained by 
numerical optimization; as for RNs the optimization 
procedure gives only slightly better results 



model devised for these systems on more general and 
abstract networks. Our analysis has focused on the 
subtle interplay between the global network structure 
and the arrangement of the two primary light harvesting 
complexes (LHl and LH2) and their impact on the 
transport efficiency rj. We have investigated well-known 
regular structures (Cayley trees and Regular Hyper- 
branched Fractals), as well as randomly- generated ones, 
and we have considered different network-theoretical 
centrality measures in order to select different LHl and 
LH2 arrangements on the networks. 
One can identify three transport regimes, depending on 
the relation between the rate of photon capture and the 
recycling time of the reaction centers (RC) which are 
placed in the center of LHl complexes: the active, the 
transition and the saturated regime, characterized by 
the increasing average number of closed RCs found by 
an incoming exciton. 

Our results clearly show that topology is crucial for effi- 
cient transport in the systems under analysis: different 
LH1/LH2 arrangements yield relative differences 5r] in 
the transport efficiency which can be as high as 37%. 
In particular, the sensitivity of the efficiency to the 
LH1/LH2 arrangement Srj is enhanced when the number 
of LHl is lower (low light illumination conditions) and 
typically in the transition region, where 50% of the RCs 
is closed, making the LH1/LH2 arrangement even more 
relevant in this region. 

Furthermore, we are able to identify guiding principles 
for optimal LH1/LH2 arrangement, given a global 
network topology. While efficient configurations always 
correspond to lower distances traveled by excitons to 
reach open RCs, the optimal criteria to achieve this 
feature strongly depend on the topological properties of 
the networks considered. 

For regular structures, efficient configurations can be ef- 
fectively devised by means of centrality measures. Global 
(local) centrality measures induce efficient arrangements 
for networks with small (large) average distance between 
the nodes. Not only are these arrangements much more 
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efficient than a random arrangement of complexes over 
the network: they are optimal, in the sense that they 
can match the efficiency obtained through a numerical 
optimization method. These results might be of use in 
the construction of artificial light harvesting networks 
where global control over the network structure and the 
arrangement of complexes is possible. 
As for randomly-generated structures the arrangement 
which ensures the best efficiency is given by a local 
criterion (degree). However, it is not possible to devise 
arrangements which are significantly more efficient than 
the random one, neither through centrality measures 
nor trough optimization. This fact suggests that the 
transport process modeling purple bacteria membranes 
is robust with respect to limited amount of control in 
the structure design, and it is optimized to take into 
account the disorder inherent in the construction of 
randomly-assembled structures. 

Finally, real biological networks have in general equal or 
better performance with respect to the other structures 
considered. While artificial structures with optimal 
arrangement of LH1/LH2 can ensure a comparable 
efficiency in the active region, a significant gap appears 
in the transition region where real biological networks 
substantially benefit from their higher level of connec- 
tivity. 

The framework we have developed in our work can 
be extended by comprising different centrality or anal- 
ogous measures, and it can be applied to the study of 
more general processes of transport on complex networks 
with the presence of dissipation, trapping and congestion 
with the goal of deriving general design principles for 
artificial light harvesting networks. 





(a) DEG distribution CT (b) SPB distribution CT 
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(c) DEG distribution RHF (d) SPB distribution RHP 

FIG. 14: Distributions of DEG, SPB for CTs, d = 4, 
n = 4 and for RHF, f ^5,n^3 
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Appendix A: Characterization of the selected 
networks 



We report in the figures below (Fig. I14I15I16]) the 
distribution of the network-theoretical measures (DEG, 
SPB) for the selected networks. Notice that the distri- 
butions shown for random RNs (Fig. I16p are an average 
over M — 100 networks with different topology. 

As for the CTs: the DEG distribution is trivial, there are 
d{d-l)"-^ = 108 nodes with DEG= 1 (peripheral nodes) 
while all the others have DEG= d = 4; SPB is a strongly 
decreasing function of the generation of the nodes (their 
distance from the center) . The number of nodes of gener- 
ation k scales as ~ d*^, so most nodes (144) have a small 



(c) DEG distribution real (d) SPB distribution real 
LLI LLI 

FIG. 15: Distributions of DEG, SPB for real HLI and 
LLI networks 



SPB (< 0.1), while 12 nodes have SPB - 0.2, 4 nodes 
have SPB ^ 0.4, and one single node (the central one) 
has SPB - 0.8). The DEG and the SPB arc poorly cor- 
related (linear correlation coefficient c = 0.57). 
As for the RHF: the DEG distribution is simple, a num- 
ber (/ -I- 1)"~^ = 36 of the nodes (the central nodes of 
the component star graphs) have DEG= / = 5, a num- 
ber 



— 26 have DEG=: 2 and the remaining 
ones have DEG= 1; SPB is again a decreasing function of 
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(a) DEG distribution (b) SPB distribution 

m = 193 m = 193 
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(c) DEG distribution (d) SPB distribution 

m = 321 m = 321 

FIG. 16: Distributions of DEG, SPB for random 
networks 



the distance from the central node, the majority of nodes 
(166) have SPB < 0.1 while a smaller fraction has values 
SPB in the range 0.2 — 0.4, and one node (the central 
one) has a SPB - 0.8). Again, the DEG and the SPB 
are poorly correlated (c = 0.47). 

As for the RNs: the DEG distribution is peaked around 
DEG- 4, DEG-- 3, corresponding to to = 197, m = 321 
and it is skew, showing a tail for small/high DEG to = 
197, to = 321; the SPB distribution is decreasing, all 
nodes have small values of SPB (quite differently than in 
the case of CT and RHF), values of SPB (- 0.2-0.3) are 
attained only for a few nodes. The DEG and the SPB 
are poorly correlated (c ~ 0.60 for to ~ 197 and c = 0.55 
for TO = 321). 

As for the real biological networks: the DEG distribu- 
tion is sharply peaked around DEG'^ 6; the SPB distri- 
bution is irregular, all nodes have small values of SPB 
(0.01 - 0.06), with most node having SPB - 0.02. The 
DEG and the SPB display a weak correlation for LLI 
networks (c = 0.70) and a stronger one for HLI networks 
(c = 0.82). 

Appendix B: Master equation approach 

Alternatively to performing random walk simuations, 
one can evaluate many figures of merit directly from a 
Master equation approach. The single exciton transfer 
and trapping can be described in terms of the random 
migration of a single exciton in a membrane with an effec- 



tive fraction of closed RC. Indeed, for any fixed tbiock the 
number of closed RC rapidly reaches its average value, so 
that each exciton sees an effective fraction of closed RCs 
corresponding to a this average value Q. This allows 
for a Master equation treatment of both active (all RCs 
open) and saturated (most RCs closed) photosynthesis. 
More precisely, the process is governed by different prob- 
ability rates and thus can be modeled by a (Markovian) 
Master Equation (ME) 

= ^ KjnnPn (Bl) 
n 

whose solution is \p{t)) ~ e^^\p(Q)). The A^-dimensional 
vector \p) is composed by the probabilities {p„} that the 
excitons be in site n. The initial probabilities are given 
by 

{ai/a \in = LHl 
cra/cr iin = LH2 
\in^RC 

(B2) 

and transfer matrix K can be written as (9|: 

Kmn = W„in " S„inC^ Win + ^n,flC°^cs) ^ (B3) 

where: The exciton transfer rates Wmn are taken as the 
inverse of the exciton transfer times given in Table HIl 
kd.ss = t-^L^ k*^^^ = t^h^,,^, kes = t-\ In the ME ap- 
proach, the effect of the recycling time of RCs can be 
taken into account by keeping Nuock RCs closed, where 
Nbiock is determined as the average number of RCs closed 
at time t in a fixed illumination condition. For each value 
of Nhiock some relevant functionals (efficiency, exciton 
lifetime, etc.: See below) are evaluated. The procedure 
is repeated a high number of times, randomly choosing 
the RC that are closed, and one finally determines the 
average of the functionals for a given Nuock [ssj . 
In the ME approach, 77 can be simply obtained from 
the inverse transfer matrix Q. Indeed we have 

?/ = u!csit)dt where uJcs{t) is the probability of an ex- 
citon causing charge separation at time t and is given 

hy ujcs = k^sY^n^RcPj^i^) = ^csE„ei?c("|e^*b(0)> 
where \n) is the probability vector corresponding to an 
exciton being in site n with certainty. Therefore rj = 

-fc-E„6flc("l^"'|p(o)) • 

The results of the random walk dynamics have been com- 
pared with the ones obtained in the Master equation ap- 
proach and the two methods show complete agreement. 
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